Resonance plasmon linewidth oscillations in spheroidal metallic nanoparticle 

embedded in a dielectric matrix 
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The kinetic approach is applied to calculate oscillations of a surface plasmon linewidth in a 
spheroidal metal nanoparticle embedded in any dielectric media. The principal attention is focused 
on the case, when the free electron path is much greater than the particle size. The linewidth of the 
plasmon resonance as a function of the particle radius, shape, dielectric constant of the surrounding 
medium, and the light frequency is studied in detail. It is found that the resonance plasmon linewidth 
oscillates with increasing both the particle size and the dielectric constant of surrounding medium. 
The main attention is paid to the electron surface-scattering contribution to the plasmon decay. All 
calculations the plasmon resonance linewidth are illustrated by the example of the Na nanoparticles 
with different radii. The results obtained in the kinetic approach are compared with the known ones 
from other models. The role of the radiative damping is discussed as well. 

PACS numbers: 78.67.Bf; 68.49.Jk; 73.23.-b; 78.67.-n; 52.25.0s; 36.40.Vz; 25.20.Dc; 52.25.0s; 78.47.+p. 
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I. INTRODUCTION 

The surface plasmon (SP) excitation (dipolar reso- 
nance) in metallic nanoparticles (MNs) are still of great 
fundamental interest -'^^'^ since its pronounced local reso- 
nances, the position, shape and intensity can be tuned 
over wide spectral range by varying the size and shape of 
the MN or by changing the surrounding medium. Once 
excited, plasmon oscillations can damp non-radiatively 
by absorption caused by electron-phonon interactions, 
and/or radiatively by the resonant scattering process. 
Due to their interaction with other internal degrees of 
freedom, the plasmon excitations can decay that results 
in a certain linewidth F [full width at half a maximum 
(FWHM)]. The linewidth of the plasmon resonances are 
an important parameter, since it contains the informa- 
tion about the character of interactions in the system, 
in particular, about the lifetime of the surface plasmon. 
Furthermore, the information on F would help to analyze 
the specificity of electronic decay in the MNs. 

The interest to the study the plasmon linewidth is 
maintained because it is connected to the local field en- 
hancement, an effect which increases the intensity of the 
incident light near the MNs surface by several orders of 
magnitude.''"^ Since a lot of devices incorporating MNs 
gains from this effect, F is treated as a main parameter 
in applications such as field concentration for nanopat- 
terning with nanowireSf^ plasmonic nanolithography,^ 
and near-field optical microscopy,^ astigmatic optical 
tweezers,^'' surface enhanced Raman scattering^ii etc. 

The damping of surface plasmon resonances in MNs 
has been well studied both experimentally^"— and 
theoretically^"—. Usually, both the bulk and the sur- 
face damping mechanisms play an important role in the 
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surface plasmon decay. In the MNs of smaller radii, the 
penetration depth of the plasmon field reduces and be- 
comes more localized near the surface. As a result, the 
bulk-induced loss processes play only a minor role and the 
electronic excitations generated by the surface potential 
dominate. 

In the bulk material, the damping constant F is related 
to the lifetimes of all electron scattering processes that 
are mainly due to inelastic electron-electron, electron- 
phonon, an electron-defect scattering. Landau damping, 
an excitation of bound electrons into the conduction 
band, emission of electrons j^i^i etc. However, if electrons 
are confined in nanoparticles with sizes below lOOA, the 
electron-surface scattering comes into playJ^ The surface 
acts as an additional scatterer because the mean free path 
(MFP) of the electrons becomes comparable to the size 
of the particles. For very small particles, the collisions of 
conduction electrons with the particle surface dominate, 
that results in the reduction of effective MFP. According 
to this model, the damping constant F depends on the 
particle radius R. The width of the spectral line has been 
often explained^ by introducing a term varying linearly 
with the inverse diameter of the particle. 

Kreibig was the first who introduced this classical pic- 
ture of the limitation of the MFP and found good agree- 
ment with experimental results.^ An early quantum me- 
chanical model developed by Kawabata and Kubo2i. also 
predicts a 1/R dependence of the plasmon bandwidth. 
However, their model do not treat the surface as a scat- 
terer for the electrons. The 1/R dependence refiects the 
importance of the ratio between the surface area S and 
the volume V . 

In Refs. [28l - l3ll . the semiclassical theory was used to 
evaluate the surface plasmon linewidth in MNs. The os- 
cillations of the surface plasmon lifetime as a function of 
nanoparticle size have been found with the use of such 
approaches. By numerical calculations, it was found that 
the nonmonotonous size-dependent structure in the line 
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shape can be attributed to the sheh effects. For smah 
alkahne-metal clusters has been shown that the angular- 
momentum-dependent electron-hole density-density cor- 
relations lead to an oscillatory size-dependence of the 
plasmon linewidth due to Landau damping. 

But not only these correlations and Landau damp- 
ing cause the linewidth oscillations. The purpose of the 
present paper is to calculate the linewidth of the surface- 
plasmon resonance when the MFP of the electrons is 
vastly larger than the particle size and their scattering 
on the particle surface plays the most important role. 
Within the framework of the kinetic approach, we found 
that the the surface plasmon linewidth exhibits the os- 
cillations as a function of the nanoparticle size, shape, 
and light frequency. We will demonstrate how they cor- 
respond to the oscillations detected earlier—"— in nu- 
merical calculations based on the time dependent local 
density approximation. 

Therefore, the specific interest is to study the behav- 
ior of the resonance linewidth with varying of such a pa- 
rameters as the particle radius, its shape, the dielectric 
constant of surrounded medium, the deviation of light 
frequency from the frequency of plasmon resonance and 
some others. The dielectric constant of the MN bulk is 
taken into account as well. 

The rest of the paper is organized as follows. The theo- 
retical background to the problem is presented in Section 
IL Section III contains the study the plasmon linewidth 
in some specific cases. Section IV is devoted to the dis- 
cussion of the obtained results. The radiative damping 
is studied in the Section V and Sections VI contains the 
summary. 



lowing empirical relation 



Is 



6.23.28.33 



^off 



(1) 



is often used, where IcS is reduced effective MFP and A 
is a phenomenological factor. But such a formula can 
be applied to the MNs of a spherical shape only in the 
case when the MFP of the electrons I is smaller than 
the particle size d. If the shape of MN differs from the 
spherical one and I ^ d, then the formula ([l} can no 
longer be used. 

We will consider the scattering of electromagnetic 
(EM) waves on metallic nanoparticle, the size of which 
d is much smaller than the wavelength of light A ^ c/cj, 
or kd ^ 1, where c is the speed of light in vacuum, uj 
is the angular frequency of the light and k refers to the 
wave number. Then the EM field around the MN can 
be considered to be homogeneous and across a particle 
as uniform, such that all the conduction electrons move 
in-phase producing only dipole-type oscillations. An ad- 
ditional limitation to the particle size is connected to 
the electron MFP. The MN size we take as being much 
smaller than I. Then collisions of the conduction elec- 
trons with the particle surface become the most impor- 
tant relaxation process. A diffuse boundary scattering is 
assumed to be a good approximation in this case. 

Let us consider the MN embedded in a homogeneous, 
transparent medium with a dielectric constant em- The 
average power absorbed by the MN from the external EM 
field in the dipolar approximation is given 
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II. THEORETICAL BACKGROUND 

The interaction of light with a MN embedded in the 
medium is studied in the framework of classical optics, 
assuming that the particle and the medium are continu- 
ous, homogeneous, and characterized by their dielectric 
function. To overcome the problem connected to the in- 
homogeneous line broadening (connected to the size and 
shape distribution within the particle ensemble), we will 
restricted ourselves only to a single MN, which directly 
yields the homogeneous linewidth. The single-particle 
scattering spectroscopy based on the localized surface 
plasmon resonance spectra of a single metal nanoparticle 
is well developed nowj -^^'^^'^^i^^ 

In the classical case of free electrons in bulk metal, the 
damping jb{= v) is due to the inelastic scattering of the 
electrons with phonons, lattice defects, or impurities {u 
refers to the electron collision frequency), which shorten 
the MFP. In this case, the relation 7b — vp/loo holds, 
where wj? is the Fermi velocity and loo is the MFP of 
conduction electrons in the bulk. Similarly, to estimate 
the surface effect in electron-surface scattering, the fol- 



where V is the particle volume, a'jj is the principal com- 
ponents of the real part of the conductivity tensor, Ej in 
are the components of the electric field inside the MN, 
which are connected to the components of the uniform 
external electric field (coordinate independent) Ej^^ by 
the relation^ 
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where Lj are the principal value of the j-th component 
of the depolarization tensor that is also known as a ge- 
ometric factor. The explicit expressions of Lj for a MN 
with a, particular shape can be found elesewhere (see, e.g., 
Refs. [ail, [3, and dl). By e(= e' + e"), the complex di- 
electric permittivity of the particle material, which would 
correspond to a given frequency w, is denoted. Substi- 
tuting Eq. © in Eq. ©, one finds 
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Here d and e" are, respectively, the real and imaginary 
parts of the dielectric permittivity of the particle mate- 
rial; ej„ and are, respectively, the real and imaginary 
parts of the dielectric permittivity of the matrix. 

To account the surface effect, the dielectric 
functioit^ii^^ 



e'(w,4ff) = einter(w) + 
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(5) 
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is used often, where eintcr accounts for the mterband elec- 
tron transitions, and the expression in the parentheses 
applied for the miraband electron transitions. The effect 
of the surface is reduced simply to the addition of a term 
7s = 7s (/off) in the denominators of Eq. ([5]) in the form 
of Eq. ©. 

In our model we will seek the expression for 7s di- 
rectly from the kinetic equation method. We restrict 
ourselves to the case of frequency close to the frequency 
of the bulk plasma oscillations of electrons in metal 
[wpi = (47rne^/m)-^/^]. Then both terms 76 and 7s can 
be neglected in the denominators of Eq. ([5]) and the di- 
electric function for a free electron gas can be expressed 
within the Drude-Sommerfeld model as^i^i^>^ 
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The effect of a core polarization (the miraband contribu- 
tion) in these formulas has been neglected as well. For 
certain metals, e.g., K, Na, Be, Al, Mg, etc., this is a good 
approximation for the frequency range of interest. But 
for other metals, e.g., Ag, the interband transitions or the 
core polarization ought to be taken into account. Then, 
the core polarization replaces the unity on the right-hand 
side in Eq. dH) (for e') by 1-1- einter- 

For the sake of simplicity, we will assume that the di- 
electric matrix has no influence on the MN and can be 
characterized by 



j(a;) — const 
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i.e., the dielectric constant of the surrounding medium 
is assumed to be frequency independent. However, it 
may happen in actual cases that the dielectric medium is 
strongly absorptive at frequencies below the ujpi . If that is 
the case, then the is strongly dependent on frequency 
and contributes to the attenuation of the oscillations. 

Using Eqs. © and ([7]), Eq. @ can be presented in the 
form 

W'M = j4E|£^f r^Ima,,(a.), (8) 



with 
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is the imaginary part of the dynamic polarizability in 
Lorenzian form The designations introduced in Eq. © 
have the following senses: 
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is the dimensionless parameter,— 
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are the frequencies of the surface plasmon resonances, 
and 
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defines the damping rate, linewidth or, correspondingly, 
the decay time of the plasmon resonance due to electron 
scattering both from the bulk and from the surfaces of 
the particle. For understanding the decay mechanism 
of the electron plasma oscillations the knowledge of the 
decay time is of central importance. 

The /-sum rule ought to be fulfilled in these processes 



w Im oijj (lo) du) = N, 
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where N is the total number of electrons in the MN. 

In the case of medium with e„, — > 1, Eqs. (|8))-([T0)) are 
reduced merely to 



^j,m — Ij 



(14) 
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So, the decay time of the plasmon resonance is the electric 
conductivity of the MN at the light frequency (optical 
conductivity) multiplied by a geometrical factor. 

To study the dependence Tj(uj), it is necessary to find 
the real part of the conductivity tensor as a function of 
frequency. There are different possibilities to calculate 
a'jj (uj) for the different frequency region. Below, we will 
apply the kinetic equations approach. Benefit of this ap- 
proach is that it permits one to study the effect of the 
particle shape on the measured physical values. Second, 
it enables us to investigate the particles whose sizes are so 
small that the particle surface start to play an important 
role. 

Using this method, we have found earlier— the general 
relation for complex conductivity tensor 
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where e and m are, respectively, the charge and mass of 
an electron, v is the electron collision frequency, Vj is the 
j-th component of the electron velocity and vp refers to 
the electron velocity at the Fermi surface. The complex 
\[' function entering in Eq. (ITTI) . has the form 



*(<z) = $(9) 
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with 
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(7=— (iy-H, (19) 



axis, respectively. In the case of MN with a spherical 
shape, — R± = R, v' = v, and the ^ function ceases 
to depend on the angle 9. 

Below, we will consider some particular cases that en- 
able us to derive the explicit analytical expressions for 



III. PARTICULAR CASES 

Let us introduce the frequency of electron oscillations 
between particle walls as 



and v' {— (^v) is a "deformed" electron velocitya^ with 
the "deformation" coefficient c^j — R/Rj- The last sum- 
mand in Eq. ([T5|) represents the oscillation part of the ^ 
function and the first one refers to its smooth part. 

Further, we will restrict ourselves to the nanoparticles 
with a spheroidal shape only. In this case, we have found 
that the components of the conductivity tensor for light 
polarized along (||) or across (_L) the rotation axis of a 
spheroidal MN are 
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Depending on sizes of MN, its shape and temperature, 
the variety of relations between frequencies , and ujsp 
can be achieved. For example, for the Na nanoparticle 



with the radius of R < 2 A, 



On the other 
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hand, with R > 126 A, the electron oscillation frequency 
becomes i/s < 1^, where j/~4.24-10^'^s~^is estimated for 
the Na at 300° K. This leads to a different expressions 
for (j{Ld), which can be used in Eqs. (|12p and ([T6| for 
calculation of the plasmon linewidth. 



(20) 

where n is the electron concentration and 9 is the angle 
between rotation axes of the spheroid and direction of an 
electron velocity. Here and below, the upper (lower) sym- 
bol in the parentheses on the left-hand side of Eq. (^0)) 
corresponds to the upper (lower) expression in the paren- 
theses on the right-hand side of this equation. The sub- 
script V = vp means that the electron velocity in the final 
expressions should be taken on the Fermi surface. The \1/ 
function in Eq. (|20p varies now with the angle 9 because 
the q for a spheroidal particle becomes dependent on the 
angle 9, namely 



=vf A. Plasmon linewidth of a spheroidal MN with an 
account for the bulk damping 

The components of the electric conductivity can be 
represented in terms of an expansion in a power series of 
v/vs- Then from Eq. (j20p within the frequency region 
io ^ Vs ^ V, one gets 

In the case of a; 3> I's , but v ^ Vsi we have 
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where and i?^ are the spheroid semiaxes di- 

rected along and across the spheroid rotation axis, 
respectively 4S The semiaxial ratio R^/R\\ is a measure 
of the shape of MN. The semiaxes are connected to the 
radius of sphere R of an equivalent volume through the 
relation R^ = R\\R\. 

The species of Eq. (pij) for q is governed by the "de- 
formed" electron velocity entering into Eq. (jl9p . which 
in the case of a spheroidal MN takes the form 



The factors /|| and I± entering in Eqs. (l24l) and ((25)) 
depend only on the spheroid axial ratio x {— i?j_/i?||), 
and have the form 
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and 



v' — vR\ 
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y'iO), (22) 
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x^ - 1 
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4(x2 - 1) 
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(27) 



where W|| = vcos9 and fj 



; sin 9 are the velocity 



components along and across to the spheroid rotation 



For nanoparticles with a spherical shape (x — > 1), /|| = 
I± = 1/3. In the case of the particles with a pro- 
late shape (x < 1), one should perform in Eqs. (pS)) 
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and ([27|) the following replacement: In |a; + -y/x^ — l| 
i arcsinVl — x^- For strongly oblate (x ^ 1) or prolate 
(x < 1) MNs, Eqs. ([Ml) and (H?]) reduce to the form 

T ^^—(^\ T ^^^(^\ 
\l),oh- 4 \ii2)^ (l)^P™i- 16x1/3 1^3/2^' 

(28) 

respectively. 

Therefore, the line- width of the plasmon resonances for 
two light polarizations (|| and _L to the spheroid rotation 
axis) takes the form 

, , / 9 \ 
Tf II ^ w = — — — V + -vj( in H , 

(29) 

provided that uj ^ ^ and 

^pi-^(i) 
^(")^''^" <n + i(i)(l-e;„) 

X (,^;^Ti:;2-r- (.2 + ^2)2 ^(1) + (3o) 

I 



provided that uj ^ Vg and z/^ ^ i/. The factors L\\ and 
i_L in Eqs. ([^ and ([5(11) are the longitudinal and trans- 
verse components of the geometrical factor, respectively. 



B. Plasmon linewidth of a spheroidal MN when 
the bulk damping is neglected 

1. Plasmon linewidth of a spheroidal MN m HF limit 

The components of conductivity tensor for a spheroidal 
MN in the highfrequency limit {lo Vg) and Vg ^ can 
be represented as^ 

where R±{= Rx^^^) is spheroid semiaxis directed across 
to the spheroid rotation axis, and ?/(ep) and p{ep) are 
smooth functions dependent only on the spheroid ec- 
centricity Bp — \J\ — (a prolate spheroid), or ep = 
— 1 (an oblate one): 



--rr (l — 2e^) ./l — + -r^- arcsine„, for a prolate spheroid 

^M=<! . : f V^A . ' (32) 



(1 2e2) yrr^ - ^ In (e, + yrr^) , 



for an oblate one 



+ 2ep)^l — Gp — (l — 4ep) arcsine^, for a prolate spheroid 

^ ' ^ - (1 - 2el)^l + el +3^3(1 + 4e^) In (^Cp + ^1 + e^J , for an oblate one 
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Eq. ([T2I) with accounting for Eq. ([Bl]) . takes the form 
^(•11^^) = TT Ti T-'^Cin c ^ ■ (34) 



Using Eq. (fTTj) for the plasmon resonance frequencies 
(w = Wj^sp, j = II : -L), Eq. (IMl) can be reduced to the form 

r ^ ^ 9 /^?7(x)^ 

where the product Rx^^^ represent the Rjl- In this case 
the linewidth depends on the shape of MN solely through 
the functions ?7(x) and p{x). 

One can use the asymptotic expressions for functions 
r]{x) and p{x) in the cases of both the extremely small or 



the large axial ratio: 

r vr/S -I- SttxVIG, x< 1 
r/(x) ~ <^ 

[ x/2 4- l/(4x), x>l 

(36) 

[ 37r/16 + 7rxV32, x < 1 

p(x) <^ 

[ x/4-f- (-l-F41n2x)/(8x), x>l 

If one consider the nanowires and nanorods, which can 
be reasonable approximated as prolate spheroids, then 
one can put 77 ~ tt/S and p ~ Btt/IG with a sufficient 
degree of accuracy. The depolarization coefficients in this 
case (x ^ 1) look as 

L||(x) ~ x^ [ln(2/x-x/4) - 1] , ^^(x) = [l-i||(x)]/2. 

(37) 

In the case of MNs with a spherical shape 77 = p = 2/3. 
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2. Plasmon linewidth of a spherical MN 



In the case of MN with a spherical shape, one can put 



the depolarization factor equal to L 



1/3 in 



Eqs. (1^ . ([50]) . and In the most common case, 

Eq. then reduces to the form 



(38) 



To calculate the rsph(w), one need only the function 
crgpjj(w). Let us choose for illustration the case v <^ Vg. 
The kinetic approach in this case gives the expression^ 

3 w^i 2vs . ^ ^ '^vl A 
— J^s— T 1 sm 1 ^ 1-cos — 

■(39) 



^sph 



Substituting Eq. ((391) into Eq. ([38]), we obtain 



r(a;) 



ARK uj J 2e;„ 
1 sm — 



- 1 



1 — COS — 



(40) 



Taking into account only the first term in (1401) , we recover 
at e;„ = 1 the well-knowiiSi2^i^i2a l/R dependence of F 



(41) 



As seen from Eqs. pO)) and (|4ip . the lifetime (which can 
be estimated as l/r = F) of an excitation in the MN 
depends not only on the nanoparticle radius, but also 
on the frequency (at which a given excitement is reason- 
able). For frequency that corresponds to the excitation 
of a surface plasmon in MN in a vacuum, uj = Wpi/v^, 
the following relation can be obtained from Eq. (pij) in 
energy units: 



(42) 



The same result can be derived for spherical MNs using 
the relation connecting a plasmon linewidth directly with 
a dielectric function^'"*^ 



2e"(w) 



(43) 



\de'/du}\ 



To make sure in that, it is enough to rewrite Eq. p3 
using Eq. (j6]) in the form 



An 



Wpi, 



(44) 
2/3 into 



Then substituting Eqs. (1311) with rj^ ~ pn 
Eq. (|13]), we come to the result of Eq. 

The oscillating terms in Eq. (HUl) give rise to the oscil- 
lation of F around of Fq as a function of both the particle 



radius and the frequency. They can be represented at the 
frequency of a surface plasmon as follows: 



3^3 h 



sm 



4 



R 



t^pi 
\/3 vp 



2e'„ 
1 — cos 



1 

2i?Wpi 



.(45) 



The amplitude and period of oscillations can be evaluated 
by means of the following relations 

\/3 TT-UF 



Aujp,R^{2e'^^ 



1) 



T 



LOplR 



(46) 



respectively. 



IV. DISCUSSIONS OF RESULTS 

A. Plasmon linewidth in a general case 

In order to study the significance of the oscillatory be- 
havior in more general situations, it is necessary to per- 
form the numerical calculations in Eq. (|12l) with the use 
of a general expression for the conductivity tensor given 
by Eq. (I^Hl) . By analogy with Eqs. (|40p . the expression 
([T^ with ajj{uj), presented by Eq. (pOI) . can be separated 
as well in two terms describing the smooth and oscilla- 
tory parts of the plasmon linewidth: Fq-I-Fosc- To get the 
smooth part, we restrict ourselves only to the first term 
in Eq. (|18p . To have oscillatory part, we retain only the 
last term in Eq. ([ll 




FIG. 1. (Color online) Linewidth of the surface plasmon resonance 
vs radius (in units of Bohr radius as — 0.53 A) of the spherical Na 
nanoparticles embedded in medium with = 1. The hollow rings 
correspond to the TDLDA calculation data taken from Ref. [29l . 
The smooth term ro(_R) is given by the dashed line and the solid 
line corresponds to the sum of both terms in Eq. II18I I. 

Fig. 1 shows the full linewidth F = Fq + Fqsc of the 
plasmon resonance as a function of a particle radius. 
Here and below, the results were obtained by numeri- 
cal evaluating of Eqs. ((T^ and ([^0)1 for the Na nanopar- 
ticle with the parameters'^: Ue — 2.65 x 10^^ cm~'^. 
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VF = 1.07 X 10^ cm/s and Wpi = 9.18 x lO^^ 8"^ We 
have taken i^soo^c — 4.2 x 10^'^ s~^, which, is assumed to 
be an appropriate value for the Na nanoparticle. Know- 
ing only rie, the parameters Wpi and vp can be estimated 
therewith by means of the formulas: 

27Th ( 3 \ ^^"^ 

Wpi = ^/I^m^e^Jm, vp = ( rie— I . (47) 

m \ Stt / 

Mostly, for the kinetics, the following inequality could 
be met: 

2TTh , , 

-R>i?dB = . 48 

mvp 

For the Na particles, e.g., R ^ i?dB — 7A. Since we can 
not apply the kinetic method to the range of R, where 
the quantum effects (such as, e.g., the Landau damping) 
play an important role, we have restricted ourselves to 
some minimal value of i?min, from which the R should be 
measured. To fit our calculations to the TDLDA calcu- 
lations for the Na cluster presented in Ref. H^, we choose 

though, i?min = -RdB/2. 

As can be seen in Fig. 1, the calculated smooth compo- 
nent of the SP linewidth ro(-R) is inversely proportional 
to the radius of the nanoparticle. The resonance peaks 
are sharper for larger particle radii and tend to broaden 
for lower particle radiuses. The oscillating terms rep- 
resent an important correction to ro(i?), especially at 
small particle radii. Additionally, one can see that the 
reduction of the MN radius up to R/as — 12 leads to 
the increase of F to the maximal value around of 0.78 eV 
(in the time domain it corresponds to the minimal SP 
dephasing time of T,„i„^sp ~ 1.7 fs). In other words, the 
dephasing rate of spherical MHs increases for larger par- 
ticles. A further reduction of i? to i?/aB ~ 8 causes the 
decrease of F to values as low as 0.3 eV. Such a behavior 
of F(i?) (without oscillations) resembles qualitatively the 
measured one for the Ag nanoparticlesj^ 

Our result for the Na nanoparticles mainly agrees with 
a similar results obtained in Refs. [29lf30l Experiments on 
alkaline clusters with a diameter in the range of 10 — 50 A 
in vacuum^^ yield a linewidth of the order of F 1 eV. 
Our calculated value is slightly smaller, but of the same 
order of magnitude as the experimental one. 

B. Particle shape effect 

The next interesting question is the effect of a particle 
shape on the plasmon relaxation dynamics. The shape 
factor of the ellipsoidal particles can be parameterized 
by the R±/R\\ = x ratio. Oblate ellipsoids correspond 
to the case of x > 1 (pancake shaped), spherical parti- 
cles - to .T = 1, and prolate ellipsoids - to x < 1 (cigar 
shaped). As already was outlined,— i^i^ any change of 
the nanoparticle shape from a sphere and thus introduc- 
ing of an anisotropy, results in the splitting of the sur- 
face plasmon resonance into two modes: a transverse one 



(perpendicular to the spheroid axis of revolution) and a 
longitudinal one (parallel to this axis). 

In the investigated size regime, the frequency of the 
resonance absorption is rigidly determined by the shape 
of MN*£iil With increasing particle prolateness, the ||- 
component of the surface plasmon resonance peak (at 
u! < ojsp, recall that Wgp = (-^Jpi/VS refers to the SP fre- 
quency for MN with a spherical shape) shifts to the red 
side, whereas the _L-component (at uj > ujsp) of the plas- 
mon peak shifts to the blue side of the spectrum.^S, For 
more and more oWate-shaped MNs, on the contrary, the 
_L-component of the peak (at lj < ujgp) shifts to the red 
side, and the ||-component (at uj > ujsp) shifts to the blue 
side of the spectrum. 
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FIG. 2. (Color online) Linewidth of the components Fji (thin line) 
and Fx (thick line) of the surface plasmon resonances (^||,p; and 
wx.pi vs axial ratio _Rx IK\ \ for the Na nanoparticles with the radius 
20A embedded in the medium with era = 1- The short-dashed 
(Fg ||) and long-dashed (Fo,x) hues correspond to the calculations 
of Eq. I|18|l without account for the oscillations terms. 
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FIG. 3. (Color online) The same as in Fig. 2, for the prolate MNs. 
Inset shows the F(x) at small x ratio. 

We investigate also the behavior of the plasmon reso- 
nance linewidth of a MN with the variation of MN shape. 
In Figs. 2 and 3 the dependence of the two components 
of the linewidth on the axis ratio x of the oblate and pro- 
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late Na particles, respectively, are presented for two SP 
resonance frequencies W|| and u}±^sp- The radius of the 
sphere of an equivalent particle volume was taken to be 
R = 2QA. It is supposed that the Na particle is placed 
in vacuum. 

The numerical calculations were carried out using 
Eqs. ([T2|) and (|20|) with the same numerical parame- 
ters as was above-mentioned. The smooth parts of the 
damping ro.y (short-dashed lines) and To,± (long-dashed 
lines) correspond to the case when the oscillatory terms 
in Eq. was neglected. The inset in Fig. 3 shows the 
behavior of F at small axes ratio. One can see that the 
smooth part of the damping Fqji has the maximum too. 

As is seen from Figs. 2 and 3, both components of the 
linewidth oscillate around the dashed lines with the axial 
ratio variation. The period of these oscillations depends 
on the particle volume and enhances with the x growth. 
The amplitude of oscillations enhances also with volume 
contraction. However, the oscillations of F are disap- 
peared for some axis ratio x, when the F reaches maxi- 
mum (Fig. 2). 

The peaks of F mean that among a variety of shapes 
of MN there exist such, for which the surface plasmon 
lifetime is minimal. Vice versa, the minimums of F means 
that among a variety of shapes of MN there exists such, 
for which the surface plasmon lifetime is maximal. 

The maximal and the minimal values of two linewidth 
components, deduced from Figs. 2 and 3, are tabulated in 
Table I for both the oblate and the prolate Na nanopar- 
ticle with two different radii. From Table I we observe 
that with bulking of an oblate MN (equivalent to an in- 
crease its radius) , the maximum of both components of F 
shifts to the side of the larger values of the R±/R^^; the 
absolute value of the F^-component tends to diminish 
and the F||-component retains approximately the same. 
The minimal linewidth Fmin was detected only for the 
transverse component of the F. With an increase of R, 
the location of the Fmin does not change, but its intensity 
is substantially reduced. 

In the case of prolate Na particles with bulking of a 
MN, the maximum of both components of F, on the con- 
trary, shifts toward the side of smaller values of the ratio 
R±/Rl\; the absolute value of F is reduced for the F||- 
components, but for the Fj^-components retains roughly 
the same. The minimal linewidth Fmin, was detected also 
only for the transverse component of the F. The mag- 
nitude of the Finin does not change with growing of a 
particle size, but its location is drastically shifted toward 
the side of small x. 

In general, the resonance plasmon damping in the 
oblate Na nanoparticle was found stronger along the 
spheroid revolution axis than the one across this axis. For 
the prolate Na nanoparticle, on the contrary, the damp- 
ing along the revolution axis was weaker than the one 
across this axis. This result holds regardless of whether 
the photo-excitation is close to the surface plasmon res- 
onance or far from it. 

If the MN is embedded in the dielectric media with 



Em > 1, then the environment effect ought to be taken 
into account. 



C. Environment effect 

Because the effect of an electric field on the embed- 
ded nanoparticles becomes weaker in a dielectric media 
proportionally to its refractive index, the environment ef- 
fect plays, additionally, an important role. The possibil- 
ity to measure the optical effects in any medium plays 
an important role for numerous applications in which 
the nanoparticle plasmon is used as an optical sensor 
{dr/dem) for its dielectric surrounding Persson have 
studied this problem theoretically for small silver parti- 
cles embedded in various matrixes.'^^ The spectral pecu- 
liarities of an environment effect recently were investi- 
gated for the Ag and Au nanoparticles, for instance, in 
works [i^ - lioj . Below, we will present the result of our 
calculation of the linewidth for the spherical Na nanopar- 
ticle embedded in a different dielectric media. 
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FIG. 4. (Color online) The total (solid lines) and smooth (dashed 
lines) linewidth of the surface plasmon resonance vs dielectric con- 
stant of an environment media for the Na nanoparticles with the 
radii of 20, 30, and 50A. 

The resonance plasmon linewidth, is plotted in Fig- 
ure 4 against the dielectric constant of various dielectric 
environments for the Na nanoparticles with three differ- 
ent radii. The numerical calculations were performed for 
a plasmon resonance frequency uj = Wpi/ ■y/2e^ -I- 1, with 
the use of Eqs. (fT^ and ([20)1 . As one can see in Fig. 4, the 
total linewidth of the surface plasmon resonance slightly 
rises and oscillates round its smooth part when the di- 
electric constant of the environment is increased. The 
oscillations are well pronounced for the Na nanoparticles 
with the small radii and disappeared for NP with a larger 
radii. In experiments for MNs with small sizes (where 
the bulk contribution to the F is decreased),— only the 
linewidth reduction with was fixed. 

Another results are as follows. With an increase of the 
particle radius, the linewidth of SPR appreciably falls 
and oscillates around an constant value in higher dielec- 
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TABLE I: The extremal values of the linewidth components (Fy and r±) for the Na nanoparticles with different radii. 
Oblate Na (a; > 1) | Prolate Na {x < 1) 



R{k) r. 


,max(eV) 


X 


rx,,^ax(eV) 


X 


ri,„.i„(eV) 


X 


r||,max(eV) 


X 


rx,max(eV) 


X ^^ 


_,min(eV) 


X 


20 


3.67 


102 


0.555 


15 


0.2185 


1.5 


0.375 


0.088 


2.765 


0.0005 


0.67 


0.015 


20 






0.2195" 


1.28 


0.2165 


1.03 






0.684" 


0.023 






100 


3.69 


823 


0.31 


51 


0.076 


1.5 


0.14 


0.028 


2.765 


0.8 X 10"^ 


0.67 


0.00023 


100 


















0.684" 


0.00035 







"This is for the second maximum of Fx . 



trie constant environment. The magnitude of these oscil- 
lations is the greater the smaller particle is and enhances 
markedly with e^. 

The resonance peak position gets shifted with chang- 
ing the dielectric function of the surrounding medium. 
The spectral direction of a shift depends on numbers of 
factors, well studied in earlier publications 

D. Frequency dependence of F 

Let us study how the linewidth of a surface plasmon 
resonance changes with the deviating of an incident light 
frequency from the resonance one. 

Fig. 5, depicts the behavior of a plasmon resonance 
linewidth vs the deviation of a light frequency from plas- 
mon resonance one for the prolate Na nanoparticle with 
Ri_/ R\\ = 1/10. Because there is a direct connection be- 
tween the plasmon resonance frequency ujsp and the axes 
ratio R_\_/R\\, the dependence of T on the frequency can 
be treated as well as the dependence of T on the shape 
of spheroid. 

The numerical calculation were carried out, as was 
done above, using Eqs. p2|) . (|20|) and the same numer- 
ical parameters. The frequency scale is normalized to a 
plasmon resonance frequency ujsp = Wpi/v^ = J7 for a 
spherical MN embedded in vacuum. 

Most remarkable (see Fig. 5) is the strong damping of 
plasmon oscillations at a frequencies uj much lower than 
the Wsp- This result is due only to the surface electron 
scattering. Our calculation of the plasmon resonance 
linewidth vs resonance energy qualitatively correlate with 
the its behavior measured for single Au nanospheresJ^ 
But a similar calculation for the Ag nanoparticles with 
account for the influence of the bulk dielectric proper- 
ties of the NPs on the damping process^! gives, on the 
contrary, the increase of F with the energy of a plasmon 
resonance. This may be due to the fact that the bulk 
damping exceeds the surface one for MNs with the great 
radii. 

As can be seen from Fig. 5, the _L-component of the F 
exceeds the || -component one within a broad frequency 
range. The inset shows that the damping is considerably 
decreased as the radius of the spherical particle increases. 
Additionally, the oscillations of the F gradually disappear 
as a particle becomes larger. 
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FIG. 5. (Color online) Linewidth of the components F|| (thin line) 
and Fx (thick Une) of the surface plasmon resonances for prolate Na 
particles with x = 1/10, embedded in the medium with trn = 1, as 
a function of frequency in the vicinity of the cjsp . The short-dashed 
(Fq II ) and long-dashed (Fq.x) lines correspond to the calculations 
with accounting for only the smooth term in Eq. I II8I I. Inset shows 
the same dependence for the spherical Na particles with R = 20A 
and lOOA. 

Previously, we have studied^ how the ratio between 
the transverse F_l and longitudinal F|| components of the 
plasmon resonance half-width depends on the degree of 
ellipsoid's oblateness or prolateness for frequency region 
that is located higher or lower than the characteristic fre- 
quency of an electron reflections between particle walls. 

Up to now, we considered only nonradiative processes, 
when the electron scattering in the MN dissipates oscil- 
lation energy into heat. Below we dwell shortly on the 
radiative processes. 



V. RADIATIVE DAMPING 

The problem of a damping of the electron energy due 
to the radiation of a portion of the collective electron 
oscillation energy into the optical far field has been ex- 
tensively studied in the literature/^ii^i^"~ 

As we have seen from the above sections, the plasmon 
linewidth for small MNs (i? < lOOA) essentially depends 
on electron collisions with a particle surfaces (dissipa- 
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tive damping). However, for large MNs {R ^ lOOA), 
our mechanism do not account for the dissipation of the 
electron energy due to inverse transformation of localized 
plasmons into propagating electromagnetic radiation (ra- 
diative damping). 

The relative contributions from radiative damping 
through the resonant scattering and absorption are also 
strongly depend on the particle size. In particular, it is 
known^ that the plasmon absorption is the only process 
in small particles, whereas both the absorption and the 
scattering are present in large particles, with the latter 
becoming more dominant as the particle size increases. 
The phenomenon is based on interplay of the usual dissi- 
pative and radiative damping, where the latter is related 
to inverse transformation of localized resonant plasmons 
into scattered light. 

The charges, under approach of classical electrodynam- 
ics, radiate when they move with acceleration. To cal- 
culate the line broadening that is entirely caused by an 
increase of Trad due to the radiant effect, we will use 
the time dependence of a classical dipole oscillator. The 
force of a decelerative radiation of a dipole under an inner 
electric field (see Eq. (jS])) can be presented as 



F,.ad(t) 



2e 
3? 



[l + L(e/e„,-l)]-i 



(49) 



In the case of a medium with £„ = £ = !, Eq. p9)) 
transforms into those well-known from the classical 
electrodynamicSfSi The linewidth due to the radiative 
damping of dipole vibrations is connected to the F by 
means of 



rad 



-Im 



m 



■ rad 



it) 



N, 



(50) 



where N is the number of free electrons in the MN. Sup- 
posing d(t) = do exp (— iwt), we obtain for j-th compo- 
nent of a radiative linewidth the following expression: 



r 



j,rad — „ 



3 mc^ 



Nlm 



(51) 



or with accounting for the dielectric matrix properties 
given by Eq. (O, Eq. (jSTj) is reduced to the form 



rj,rad(t^) - o 



3 mc^ 



1/2 



(52) 



Taking into account the expression ([6]) for e", we get 



rj,rad(^^) 



Stt e'^Lu LjCj'.Auj) 



3 mc^ 



N 



1/2 



(53) 



Different relations can be employed for a'^^ , depending 
on the frequency regime. For instance, in the highfre- 
quency limit, when Eq. (I3ip can be applied, we obtain 
for two component of the linewidth of a spheroidal MN 
embedded in a medium with Cm the following equation 



(I), rad 4toc3' 



(54) 



where the functions i]{ep) and p{ep) are given by Eqs. (p2)) 
and ([33]), respectively. In the case of a spherical particle, 
Eq. (|54|) with accounting for Eq. (fTTj) . can be rewritten 

as 



rad,i 



le 

6 mc 



Wpi^UF (1 + 2e 



U/2 



R 



1/2 



(55) 



The increase in linewidth from radiative damping is 
proportional to the surface area of the MN. The effect is 
weaker in a higher dielectric constant of environment. 

The estimations of the F for the spherical Au, Ag, Cu 
and Na particles with 2R = 200A embedded in water 
(e„i = 1.78) give: 0.942, 0.927, 1.82, and 1.13 meV, re- 
spectively. It is known from Mie scattering theory^ that 
only 1.5% of the total damping rate in 2R ~ 200 A metal- 
lic spheres is due to the radiative decay. 

In order to take into account radiative damping to- 
gether with collisions of free carriers with the MN sur- 
face, the effective collision frequency Fcff = F + Frad must 
be introduced. 



VI. SUMMARY 

We use the kinetic approach to study the plasmon reso- 
nance linewidth for the metal nonspherical nanoparticles 
embedded in any dielectric media. It enables one to cal- 
culate the linewidth in the case that the free electron path 
is much larger than the particle size and the scattering 
from the particle surfaces plays a dominate role. 

The general formula is proposed for a damping rate or 
a decay time due to electron scattering from the bulk and 
particle surfaces. By means of this formula one will have 
a possibility to evaluate the linewidth directly through 
the tensor of polarizability of the MN. 

The electron surface-scattering contribution to the 
plasmon damping in a simple case of a spherical metal 
nanoparticles is studied in detail. It is clearly shown that 
the resonance plasmon linewidth oscillates as the particle 
radius increases. The oscillating terms represent an im- 
portant correction to the linewidth, especially at a small 
particle radii. This result for the Na nanoparicles is in 
well agreement with the numerical time-dependent local 
density approximation calculations. 

With changing the MNs shape from spherical to the 
spheroidal one, the single plasmon resonance splits into 
two components: the longitudinal and a transverse one 
to the spheroid rotation axis. Both components of the 
linewidth oscillate with the axial ratio altering. The pe- 
riod of these oscillations depends on the particle volume 
and is enhanced with an increase in the axial ratio. The 
amplitude of oscillations enhances also with the contrac- 
tion of the particle volume. The behavior of the linewidth 
extrema are studied for both components of an oblate and 
a prolate MNs with different volumes. 

The resonance plasmon damping in the oblate Na 
nanoparticle was found stronger along the spheroid revo- 
lution axis than the one across this axis. For the prolate 
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Na nanoparticle, on the contrary, the damping along the 
revohition axis was weaker than the one across this axis. 

The size-dependent oscillations of the linewidth also 
depend on the dielectric constant. For the first time, we 
detect that the oscillations of the surface plasmon reso- 
nance may occur also with an increasing of the dielectric 
constant of the surrounding medium. The oscillations are 
well pronounced for nanoparticles with the small radii 
and disappeared for NP with a larger radii. The magni- 
tude of these oscillations is the greater the smaller parti- 



cle is and enhances markedly with e„i- 

The effects of both the particle shape and the environ- 
ment on the plasmon resonance linewidth are illustrated 
by the example of the Na nanoparticles with a different 
radii. 

The contribution of the radiative plasmon decay is dis- 
cussed as well. 

Our theoretical results should be important for the 

analysis of the transport and optical properties of MNs 
under an exposure of short and strong laser excitations. 
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